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A new gauge fixing that imposes the U(l) real and quaternion real condition on the correlator of 
domain wall fermion on the left wall and right wall is proposed. The method is applied to the 
calculation of the mass function, the QCD effective coupling and the charge form factor of the 
proton. By assuming topological charge change \Q\ = 1, the charge form factor of the proton is 
approximately fitted by the dipole form factor. 

The difference of the Kugo-Ojima color confinement parameter in quenched simulation and in 
unquenched simulation is expected to be due to an anomaly caused by the vertex of two fermion 
zero modes and the Bose ghost/gluon. 
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1. Introduction 

The infrared (IR) QCD is characterized by the color confinement and the chiral symmetry 
breaking. We performed ^ lattice simulations of the domain wall fermion which preserves the chi- 
ral symmetry in the zero-mass limit using the full QCD gauge configurations of the RBC/UKQCD 
collaboration of 16 3 x 32 x 16 lattice [Q] and compared with results of staggered fermion of the 
MILC collaboration [J3p . In these lattice simulations and in comparison with other works, we 
observed qualitative differences between quenched and unquenched simulations in the QCD effec- 
tive running coupling, and in the Kugo-Ojima parameter [Q] which is calculated numerically by 
the inverse of the Faddeev-Popov operator whose zero-modes subtracted at the zero momentum 
[@, ^ 0, §]. The difference of the quenched and the unquenched simulation could be attributed to 
the fermion zero mode. In finite temperature hamiltonian QCD, zero-modes are localized around 
monopoles whose position is defined by the boundary condition on the time axis [||]. 

In Sect.2, we review the IR QCD and Gribov-Zwanziger-Kugo-Ojima theory, and in Sect.3, 
we show the method of gauge fixing using the boundary condition of fermions on the domain wall. 
In Sect.4, the quark-gluon coupling is presented and in Sect.5, the charge form factor of the proton 
is presented. Discussion and comments on the Kugo-Ojima parameter are given in Sect.6. 



2. The IR QCD and Gribov-Zwanziger-Kugo-Ojima theory 

The Landau gauge QCD or the Coulomb gauge QCD in the IR region suffers from the problem 



of gauge uniqueness. Gribov [10] defined the so called Gribov region £2, which is convex and 



bounded in every direction and every gauge orbit passes through it at least once. Zwanziger [11] 
showed that the cut-off of the functional integral at the Gribov horizon dQ. may be replaced by 
the Boltzmann weight. He defined the free field action Sq[A] and an ellipsoid in A space Q [A] and 
replaced the measure inside the Gribov horizon 

dAexp(-So)8(c- Q[A]) — > d\i y = dAexp(-So)exp(-yVQ[A}) 

where y is defined from the self-consistency (<2[A]) r = c, where the l.h.s is measured with use 
of dpLy. When 7^0, there appears in addition to the Fermi ghost c whose correlator (cc) yields 
the Faddeev-Popov determinant, but the Bose ghost <p which couples to the quark and yields the 
linearing rising potential [Oj. 

Kugo and Ojima [Q] assumed y = and the BRST symmetry in the IR region and defined the 
color confinement parameter u{k 2 ) as 

(D^ a (A v x c) b ) k = (5 MV - ^L)8 ab u(k 2 ), 

and m(0) = — 1 as the color confinement condition. In this theory, the ghost dressing function F(k) 
is expressed as 

F(k) 1 



\ + u{k 2 )+o{k 2 ) 
and the ghost propagator should be infrared divergent. 
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Quenched lattice simulations indicate, however, u(0) is about -0.8 and unquenched lattice 
simulation indicates that m(0) ~ — 1 [jg], but the ghost propagator does not show the infrared diver- 
gence. Recently, violation of the BRST symmetry and infrared finiteness of the ghost propagator 
are proven [0[ 14, 16], but whether w(0) becomes -1 in unquenched simulation is not clear. 

In perturbative QCD (pQCD), the ghost propagator is usually calculated via Fermi ghost-gluon 
coupling, but when 7^0 there is a coupling of gluon to Bose ghost which couples to quarks. 
The fermionic zero mode should modify the measure d\l and we want to evaluate the zero mode 
effects using a specific gauge fixing. 




Figure 1: The ghost-gluon loop contribution to 
the ghost propagator. Solid curves are gluons 
and dashed lines are Fermi ghosts. 




Figure 2: The quark loop contribution in the 
ghost propagator. Thick solid curves are quarks, 
zig-zag lines are Bose ghosts. 



3. The gauge fixing using the boundary condition of fermions on the domain wall 

On the lattice, the expectation value of the quark propagator S(p) consists of spin dependent 
p part and spin independent S3 part. The propagator is written as 

-is*p + ® -i^p + @ 
S(p) = 3- = +. (3.1) 

efijP + JtJ?*) sip 1 + Jt& 

We specify the propagator of the left-handed quark by the suffix L and the right-handed quark by 
the suffix 7?. With use of the renormarization factor Za and Zg, they are defined as 

n(Z(p,s)P L/R V(p,s))=Z B (p)(2N c )a L/R (p,s), 

and 

Tr(x(p,s)ipP L/R y(p,s)) = Z A (p)/(2N c )ip*/ L/R (p,s), 

where pt = - sin -^Ei (n,- = 0, 1 ,2, • • • , njA). The mass is defined as 
a rij 

U > Re[M R (p,L s /2)} yP> Re[^ L (p,L s /2)} 

3.1 The U(l) real condition 

Since we adopt the representation of fermion spinors such that the 75 is diagonal, the fermion 
field y can be gauge transformed in the fifth dimension as 

y^e'^y, iff -» ye'^ , (3.2) 
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such that on the left and on the right domain wall, the propagator becomes approximately real by 
choosing [jl|] 

I i8L in i |2 i I idR -in i i2 

\e e 1 — \\ +\e e ' — 1| 

be minimum. 

We sample wise measure the propagator by using the transformation matrix from the left 
boundary 

ciq b \ 
c do } 



a° x G\ + cP (72 + a° z a 3 + il 



which is given by the sample average of color diagonal components to the right boundary 



c Ls -i d Ls - 



3.2 The quaternion real condition 



As the transformation matrix of the fermion at s = to s = L s — 1, we adopt the ansatz [15], 



80 



e- v \( C 1 P \( ^ 
e v J y £~ l ) ^ e-v 

>C> f(Y,C)\ 
e-^- 1 i ' 



(3.3) 



In our five-dimensional domain wall fermion case, y = pL — V and /J. , v contain the phase in the fifth 
direction it] . 

2jU = icOz/n 2 - if] = (p x + ip y )£, + ip t - P z ~ m (3-4) 
2v = i03\l%\ + ir\ = (pix - ip y )C, + ip t +p z + irj (3.5) 
The quaternion reality condition of the transformation matrix g(y, Q gives 



d s -i ) I C~ l e- y J I £"V I I c do 



(3.6) 



where / = /(ft__). 



The function / and / taken in [ 15 ] is. 



doet--^e-v -X-e^- 



In general Co and cl s -\ are polynomials of £, and they satisfy c^-i^ 1 = co£ . We define 

V^ = c_ 1 C" 1 +ciC 1 +5 
where c_i = cz, _i and ci = cq and 5 is a constant, which is defined later. 
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The difference 

— tefc;Ur«4,)-(rA)(:8„ » 

should be small, if the left wall and the right wall are correlated by the self-dual gauge transforma- 
tion. 

We choose = y^-j" to make (AL/R)2i consistent with and obtain e 1 and 8, which 
approximately satisfy the quaternion reality. Corresponding to the sample-wise selection of larger 
absolute value of 8, i.e. a large zero-mode component, we assign a parameter ind = 1 for larger and 
2 for smaller, and multiply the phase (— \) md to the expectation value of the quark wave function 
that is used in the calculation of the effective mass of the quark. This phase factor can be absorbed 
in the phase e ,r] introduced to approximately satisfy the U(l) real condition. 



4. The quark-gluon coupling 

In Fig.l, the quark-gluon coupling in the Coulomb gauge of the domain wall fermion is com- 
pared with that of the staggered fermion. Red points below 1.2GeV are the experimental data of the 
JLab collaboration and the dark point is our result of momenta q = 0.7GeV/c. The dash-dotted line 
is the two loop pQCD result and the dashed line is the pQCD with A 2 condensates effect derived 
from the fitting the running coupling ^ calculated by using the gauge configuration of full-QCD 
staggered fermion of the MILC collaboration [f||. 




0.4 -0.2 0.0 0.2 0.4 0.6 0.8 



Log 10 [q(GeV) ] 

Figure 3: The effective coupling a s , gi (q) of do- 
main wall fermion and the gluon without re- 
ality selection (green^jand with reality correc- 
tion(magenta) in cylinder cut region. Dark low 
momentum points are that with reality correc- 
tion. 
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Figure 4: The charge form factor of the proton. 
Dashed line is the dipole fit of 1 /(l +q 2 /M 2 ) 2 
with M 2 = 0.7 lGeV 2 . 



5. The charge form factor of a proton 

We measure the charge form factor of a proton as a ratio of the three point function and 
the two-point function in momentum space. We calculate expectation values of etc by 
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diagonalizing the matrices whose coefficients of 0\,02 and 03 are real, and then evaluate e.g. 

{{u a + Cy 5 d b _)u c + | 7 47 M I «Cy 5 ^)<) as 

G NVN c+ {p) G NV \-{p) , G NVN h -(p) 
G»?(p) + G N a N (p) + G^UP) 1 J 

The three point function G^ N (p) is a product of 

£ 5(c + ,c +/ _,/j)(y(c +/ _,c'_ /+ , j p))7 5 5(c' f c / _ /+ , j p) t 7 5 %,^, j p). 

a,b,c,d 

where S(a,b,p) is the product of the propagator of the quark a and that of b. 

S(a,b,p) = 5(a_,c?_,p)7 5 5(ft_,ft_,/7) t 75 

In the two point function G™ (p), the factor before S(a,b,p) is replaced by the average of S(c + , c + /_ , p) 
and S(c' + c'_, + ,p)^ . 

In the product of S(c + ,c + /_,p)ysS(c' + c'_, + , pyY$ we assume contribution of self-dual field 
such that one propagator is that of left-handed and the other is right-handed. It means that we con- 
sider a change of the topological charge of one unit in the propagator. Of course it is an assumption 
which should be verified by comparison with experimental data. 

The sum of the left-handed and the right-handed contributions are to be compared with the 
vector current form factor. The form factor at p = is dominated by the V(c + /_,c , _, + ,p) of 
the left-handed fermion, while at finite p, the left-handed and the right-handed fermion contribute 
equally. The normalization at p = is 1/36 corresponding to (3 colors x 2 spin projections) 2 , and 
relative normalization at p / in the figure is multiplied by 2. (The sum of left-handed and right 
handed are close to the twice of the left-handed.) 



6. Discussion and comments on the Kugo-Ojima parameter 



In the quenched calculation, whether the u(0) becomes -1 or not depends on the treatment 



of the total derivative term [16]. On the lattice, however, this term is irrelevant. We think the 



difference of [13] and [ |T6| ] exists in another kind of total derivative term. When the quark field 



contribution in the quark-quark-Bose ghost/gluon(0 /A) vertex Q12Q 



is replaced in the zero-mass limit by the product (U i 'd^U)(U^d v U) where U are the quaternions 



of the fermion zero-modes, it becomes reminiscent of the anomalous current given in []17[], which 
is expected to stabilize the IR fluctuation: 



Jk 



da 

32n 3 J Jez 



_jabe F a 
2 U 1 KVV 



IX,Z 



where e K is an arbitrary vector. The anomaly would be cancelled by that in the Baryon sector fll8|]. 

We think that the contribution of the Bose ghost that couples to quark zero-modes cannot be 
ignored in the unquenched Kugo-Ojima parameter calculation, and that the Nature of IR QCD 
possesses a kind of supersymmetry. 
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